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The  purpose  of  this  effort  was  to  investigate  the  feasibility  of  employing  a 
piezoelectric  transducer  for  seismic  communications  applications  in  mine  environments. 
The  transducer^- a  GFE  device,  incorporated  a  high  Q,  spring-mass  resonator  to  load  the 
air  side  of  the  piezoelectric  crystal.  The  resonant  frequency  was  seismic  load  de¬ 
pendent  and  fell  in  the  range  of  300  to  340  Hz.  The  effort  included  a  mathematical 
analysis  of  the  transducer,  with  corroborative  laboratory  tests,  together  with  actual 
communications  experiments  performed  first  in  an  open-formation  limestone  quarry  and 
finally  in  an  operating  iron  mine.  This  work  showed  that  the  transducer  could  be 
coupled  efficiently  to  hard,  competent  rock  media  with  power  conversion  efficiencies 
of  70  to  80  percent  (electrical  to  mechanical)  in  typical  installations.  In  the  iron 

mli.e,  when  radiating  peak  power  was  about  30  W  (acoustical)  in  a  gated  pulse  mode, 

signals  were  received  at  a  distance  of  j60  ft.  through  a  mas rive  rhyolite  porphyry 
formation.  This  represents  a  15$  radiation  efficiency  with  an  input  of  200  watts. 
Despite  the  high  ambient  noise  level  due  to  operating  machinery,  a  signal-to-noise 
ratio  of  about  20  dB  wav  achieved  using  a  very  narrow  band  filter  to  process  the  re¬ 
ceived  signal,  and  a  data  rate  of  about  wo  pulses/sec  was  demonstrated. 
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SECTION  I 


INTRODUCTION 


The  purpose  of  the  work  reported  herein  was  to  investigate  the  feasibility  of  employing  piez  > 
electric  transducers  for  seismic  communications  and  exploratory  mapping  purposes  in  mine  environ¬ 
ments.  The  emphasis  upon  mine  environments  implies  that  conditions  of  interest  are  those  involvi  ng 
seismic  transmission  paths  or  at  least  one  terminus  of  such  a  path  deep  within  the  interior  of  massive 
rock  (or  coal)  formations.  This  distinction  is  of  significance  since  seismic  propagation  by  surface 
waves  is  thus  largely  precluded. 

The  efficacy  of  a  convenient  exploratory  mapping  system  in  mine  operations,  tunnel  driving, 
etc.,  is  rather  obvious;  correspondingly,  a  means  for  communications  through  perhaps  many  hun¬ 
dreds  of  feet  of  rock  could  be  of  tremendous  benefit  in  mine  emergencies  wherein  entrapped  miners, 
otherwise  cut  off  from  normal  communications,  could  communicate  with  rescue  teams  elsewhere  in 
the  mine  and  on  the  surface.  Moreover,  such  a  means  might  be  of  general  value  in  mine  and  tunnel 
driving  operations,  permitting  the  transmission  of  operational  messages  into  the  far  reaches,  unen¬ 
cumbered  by  the  need  for  pulling  in,  hanging,  and  protecting  vulnerable  wire  lines.  Another  impor¬ 
tant  application  would  be  that  of  providing  redundant,  damage-resistant  channels  of  communication 
in  defense  applications,  such  as  underground  command  headquarters,  etc. 

Other  applications  beyond  those  mentioned  are  also  conceivable. 

The  investigation  reported  here  is  not  to  be  regarded  as  a  comprehensive  treatment  of  the  prob¬ 
lem  of  adapting  piezoelectric  transducers  to  the  purposes  mentioned  above.  Rather,  this  was  a  some¬ 
what  specific  investigation  in  which  a  government-furnished  transducer  was  employed.  The  results, 
then,  reflect  the  performance  capabilities  of  this  transducer,  which  will  be  described  later  in  detail, 
although,  through  analysis  and  evaluation,  it  appears  reasonable  to  extrapolate  these  findings  to  a 
general  class  of  piezoelectric  transducer.  The  government-furnished  transducer,  which  was  designed 
and  fabricated  by  others,  is  a  highly  resonant  or  high  Q-type  device,  operating  at  a  frequency  of  300 
to  340  Hz  but  not  lending  itself  to  the  generation  of  sharp  impulse  signals.  Consequently,  its  useful¬ 
ness  in  exploratory  mapping  work  is  very  dubious,  since,  for  such  work,  a  very  sharp  impulse  or  high 
frequency,  gated  pulse  propagated  into  the  medium  is  normally  required.  Only  in  this  way  can  the 
return  echoes  from  various  layers,  inclusions,  etc.  be  resolved  timewise  from  the  transmit  pulse  and 
presented  for  interpretation.  Continuous  wave  exploratory  mapping  techniques  are  theoretically 
possible,  of  course,  but.  because  of  the  complexity  of  multimode  propagation  phenomena,  such 
methods  are  not  considered  generally  practicable.  This  does  not  imply,  however,  that  piezoelectric 
transducers  could  not  be  employed  for  mapping  applications,  for,  in  fact,  a  suitably  designed,  broad¬ 
band  unit  capable  of  propagating  sharp  impulse  signals  would  probably  be  the  preferred  selection  for 
such  work.  Such  units  have  indeed  been  developed.  At  any  rate,  the  transducer  employed  herein 
did  not  appear  suitable  for  the  mapping  application;  therefore,  the  effort  was  centered  upon  the  com¬ 
munication  application. 

The  investigation  included  theoretical  analysis  of  the  transducer,  backed  by  corroborative  labor¬ 
atory  tests,  for  the  purpose  of  arriving  at  a  comprehensive  understanding  of  the  performance  of  the 
unit  functioning  as  a  source  into  rock  media.  Operating  characteristics  such  as  efficiency,  bandwidth, 
radiated  power,  impedance  matching  to  the  rock  interface,  etc.  were  of  primary  interest.  This  work 
was  followed  by  communications  experiments  in  both  an  open-formation  limestone  quarry  and  an 
operating  iron  m:ne.  In  this  work,  tests  were  performed  to  determine  operating  range,  signal-to-noise 
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rutios,  effect  of  fructured  formations,  potential  data  rates,  propaptlon  velocity  and  attenuations,  etc. 
These  various  characteristics,  taken  together  with  practical  considerations  icintcd  to  equipment 
design,  etc,,  constitute  the  sort  of  information  required  to  evaluate  the  feasibility  of  this  communica¬ 
tions  method, 

To  provide  some  guidance  to  this  effort,  a  brief  survey  of  the  open  literature  was  performed, 
Owcn(l>*  presents  a  rather  comprehensive  theoretical  treatment  of  seismic  propaptlon  in  layered 
media,  with  the  added  considerations  of  controlled  seismic  sources,  received-signal  processing,  com¬ 
munications  systems  concepts,  etc,,  together  with  a  rather  complete  survey  of  the  open  literature  up 
to  1964,  In  the  vast  literature  covering  all  aspect  of  seismic  propagation  and  the  design  and  appli¬ 
cation  of  piezoelectric  transducers,  there  arc  relatively  few  references  directly  pertaining  to  this 
investigation,  As  cited  by  Owen,  the  first  efforts  of  significance  employing  a  controlled  seismic 
source  were  reported  In  a  scries  of  papers  by  IWitorf2  in  the  period  1951  to  1956,  He  used  a 
broadband  clcetrodynamlc  source  which  exhibited  very  poor  coupling  efficiency  but  successfully 
transmitted  gated  pulse  signals  in  the  100-  to  600-Hz  range,  over  distances  of  several  hundred  feet, 

In  chalk  mines  and  tunnels,  and  on  stratified  surface  layers,  Pursey(*,n>  and  Miller  and  Pursey(9',<n 
provide  theoretical  formulations  for  the  field  and  radiation  impedance,  partition  of  energy,  and 
power  radiated  by  a  source  on  the  surface  of  u  semi-infinite  isotropic  solid;  thus,  they  laid  the 
groundwork  for  an  understanding  of  the  coupling  and  propaptlon  problems.  Later,  ikrath, 
et  al(l2',5>  reported  a  series  of  investiptions  in  ihe  1963  to  1966  period  employing  a  resonant-type 
electrodynamic  source  (most  work  was  done  at  HO  Hz),  with  propaptlon  along  surface  layers  of  the 
earth,  ice-covered  lakes,  etc.  The  coupling  efficiency  of  this  transducer  was  quite  good,  and  satis¬ 
factory  signals  were  received  at  distances  of  hundreds  of  meters.  Further  work  with  transmitter 
arrays  was  also  reported,  Zcitz  and  Pakiser*l6)  cover  briefly  and  incompletely  the  application  of  a 
sonar-type  transmitter  for  exploratory  mapping  work  (grease  and  water  coupled  to  the  earth  and, 
operating  at  frequencies  of  6,  1 1,5,  and  16  kHz). 

Reprding  pcizoclcctric  transducer  design,  McMastcr,  et  al(,7)  for  example,  did  extensive  work 
in  high-power  applications,  covering  such  features  us  impedance  matching,  crystal  prcloading,  cool¬ 
ing,  etc.  for  units  operating  up  to  peak  power  outputs  of  30  kW  at  frequencies  of  10  kHz  or  higher. 

It  should  be  noted,  however,  that  none  of  the.*:  previous  efforts  has  been  concerned  with  the  design 
and  employment  of  piezoelectric  transducers  for  the  specific  application  of  communication  through 
rock  media.  However,  it  is  also  evident  that  tliis  application  is  one  of  adaptation  only,  since  the 
various  elements  of  seismic  propagation  and  piezoelectric  transducer  design  have  reached  a  high  level 
of  technical  maturity  over  the  past  years.  Accordingly,  the  chief  concern  of  the  effort  reported  here 
was  to  shed  light,  both  qualitatively  and  quantitatively,  on  the  feasibility  of  Just  such  an  adaptation. 


♦Supericript  number*  In  parentheses  re  fir  to  Htfnitnm  «t  the  end  of  thli  report, 
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SECTION  II 


PERFORMANCE  ANALYSIS  OF  TRANSDUCER 


1.  INITIAL  ANALYSIS  AND  EQUIVALENT  CIRCUIT 

So  that  a  comprehensive  understanding  of  the  behavior  of  the  piezoelectric  transducer  could 
be  acquired,  a  mathematical  analysis  of  this  unit  was  performed.  The  initial  analysis  effort  was  aug¬ 
mented  by  results  of  corroborative  laboratory  tests  on  the  transducer  so  that  a  quantitative  mathe¬ 
matical  model  was  constructed.  From  this  model,  it  was  then  possible  to  calculate  performance 
characteristics  such  as  radiated  power,  bandwidth,  and  efficiency  as  a  function  of  arbitrary  seismic 
load.  These  results  in  turn  provided  the  basis  for  the  rational  application  of  the  transducer  in  com¬ 
munications  experiments  in  actual  field  environments.  A  discussion  of  the  development  of  the 
mathematical  model  follows  after  a  brief  description  of  the  physical  characteristics  of  the  transducer. 

The  piezoelectric  transducer,  two  of  which  were  furnished  by  the  government  for  this  work, 
was  manufactured  by  the  Redford  Corp.,  Latham,  N.Y..  under  the  designation  of  Model  2427,  and 
it  has  been  further  designated  as  RADC  Model  T-100.  In  this  report,  the  unit  will  be  referred  to 
simply  as  the  RADC  transducer.  Figure  1  depicts  an  idealized  cross  section  of  the  transducer.  The 
piezoelectric  crystal  stack  consists  of  six  Vemitron  56131  (PZT-4  ceramic)  crystals  bonded  together 


Mis 

(1251b  Steel  Cyl.) 


Protective  Enve|ope 


Main  Spring 
(Al.  Alloy  Tube) 


■  I 

*  Seismic  Medium 


FIGURE  1 .  RADC  TRANSDUCER,  IDEALIZED  CROSS  SECTION 
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and  connected  in  electrical  parallel.  They  are  polarized  and  driven  in  the  thickness  mode,  which  is 
parallel  to  the  long  axis  of  the  transducer.  One  end  of  the  crystal  stack  is  held  firmly  to  the  rock 
interface  via  the  drive  plate,  using  anchor  bolts  or  other  appropriate  means,  while  the  other  end  is 
loaded  by  the  spring-mass  system.  When  the  crystal  is  driven  at  the  resonant  frequency  of  the  spring- 
mass  system,  the  mechanical  impedance  at  the  spring-mass  crystal  face  becomes  very  high  (provided 
the  Q  of  the  resonance  is  high),  thus  permitting  more  of  the  available  energy  to  propagate  into  the 
rock.  The  spring-mass  system  can  also  be  viewed,  heuristically  at  least,  as  a  quarter  wavelength  trans¬ 
mission  line.  It  is  this  resonant  drive  feature  that  permits  the  unit  to  couple  efficiently  to  the  seismic 
medium.  The  resonant  frequency  is  in  the  range  of  about  300  to  340  Hz,  a  value  which  is  sensitive 
to  the  installation  and  which  cannot  be  altered.  This  resonant  frequency  range  was  probably  selected 
by  the  designer  ■’s  a  compromise  between  radiated  power  capability  and  anticipated  signal  attenua¬ 
tion  rates  in  seismic  media;  results  of  the  effort  heie  appear  to  support  that  selection  generallv.  It  is 
mandatory  that  a  piezoelectric  crysta1  never  be  allowed  to  operate  in  tension;  thus,  a  compressive 
preload  must  be  applied  that  always  exceeds  the  maximum  operating  peak  tension  load.  This  preload 
is  provided  by  the  rod  shown  in  Figure  1,  and,  for  this  design,  this  same  loaded  rod  is  used  to  couple 
the  crystal  stack,  spring,  and  mass.  Since  the  peak  dynamic  load  of  the  unit  is  on  the  order  of 
5000  lb,  the  rod  preload  must  exceed  this  value  by  a  comfortable  margin.  Similarly,  the  transducer 
must  be  affixed  to  the  seismic  medium  with  a  load  exceeding  5000  lb  to  assure  continuous  contact. 
According  to  the  manufacturer’s  manual,  terminal  voltage  on  the  crystal  is  limited  to  1500  V  rms. 

In  envelope  dimensions,  the  unit  is  about  24  in.  long  and  7  in.  in  diameier,  weighing  about  150  lb. 

From  the  physical  configuration  of  the  transducer  described  above,  a  functional,  mechanical 
schematic  of  the  unit  is  suggested,  as  shown  in  Figure  2  along  with  a  definition  of  symbols.  It  is 
noted  that  the  crystal  stack  is  depicted  simply  as  a  force  generator  in  which  the  force,  Fc,  is  applied 
to  the  external  mechanical  elements.  (The  spring  rate  of  the  protective  envelope  is  negligible  com¬ 
pared  to  other  springs  in  the  system.)  By  examination,  it  is  possible  to  write  the  following 
expressions: 

Fc  =  wit  +  (x  +  y2  )br  +  (a'  +  y2  )kr 

Fc  =  0’i  ~x)b  +  O,  -x)k  (1) 

Fc  =  y2bs  +  y2  ks  +  (x  +  y2  )br  +  (x  +  y2  )kr 

where  overdotting  of  a  quantity  indicates  a  time  derivative.  It  is  convenient  to  work  with  the  con¬ 
cept  of  mechanical  impedance  defined  as  the  vector  ratio  of  force  to  velocity.  Accordingly,  assum¬ 
ing  sinusoidal  time  variation  with  radian  frequency,  to,  and  defining 

Vi  =y  i  =/cjy 

>'2  ~  y  2  =/«>  2 

Vm  =  x  =  jux 

Eqs.  (1 )  can  be  rewritten: 
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(2) 


rr.  ~mass 

k  ~  spring  rate  constant,  mainspring 
b  ~  viscous  damping  constant,  mainspring 
kr  ~  spring  rate  constant,  preload  rod 
~  viscous  damping  constant,  preload  rod 
ks  ~  equivalent  spring  rate  constant,  seismic  medium 
bs  ~  equivalent  viscous  damping, 'seismic  medium 
x  ~  displacement  of  mass  ; 

y  j  ~  displacement  of  upper  crystal  face 
y2  ~  displacement  of  lower  prystal  face 
Fc  ~  crystal  output  force 

FIGURE  2.  TRANSDUCER  MECHANICAL  SCHEMATIC 

we  further  postulate  an  electrical  to  mechanical  analogy,  such  that 

voltage  — *■  force 
cunent  — *•  velocity 
charge  — *■  displacement 
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5W1PB 


which  also  implies,  derivatively, 


inductance  - 
capacitance  ■ 
resistance  - 


*  inverse  spring  rate 
viscous  damping 


Then  we  can  transform  the  mechanical  re  preset  ration  into  a  sort  of  an  equivalent  “electrical”  circuit 
as  follows: 


vs  +  V2 


-U  ib  ^ 


•'I  ~vm 


In  the  Appendix,  the  equivalent  circuit  development  for  the  piezoelectric  crystal  stack  is  presented, 
and  it  is  quite  easy  to  incorporate  that  result  directly  into  the  above  development.  The  resulting 
overall  equivalent  circuit  for  the  transducer  operating  into  a  seismic  medium  !'s  shown  in  Figure  3.  It 
should  be  carefully  observed  that  all  quantities  to  the  left  of  the  ideal  electrical-to-mcchanical  trans¬ 
former  are  electrical,  while  all  quantities  to  the  right  are  mechanical.  No  loss  in  meaning  is  to  be 
implied  by  the  substitution  of  R,  L,  and  C  designators  to  the  mechanical  elements  on  the  right  side 
in  place  of  6,  m,  and  ljk  designators,  respectively.  Note  also  that  a  parallel  representation  of  the 
seismic  load  is  employed.  A  simple,  lumped  constant  representation  for  the  seismic  load  is  valid  in 
general  only  for  a  very'  narrow  frequency  range,  in  this  case  in  the  vicinity  of  resonance.  Tbs  parallel 
resistance  element,  Rp ,  is  to  be  regarded  as  the  radiation  resistance  of  the  seismic  medium,  but  this 
also  accounts  for  acoustic  energy  radiated  from  the  resonant  mass  element  and  frictional  losses  that 
might  occur  in  the  source-medium  interface  region  as  a  result  of  nonlinear  coupling  phenomenon. 
For  coupling  to  “hard”  rock  media,  'ids  latter  is  considered  negligible  so  long  as  stress  levels  in  the 
rock  do  not  become  excessive.  Net.:  also  that  the  reactive  part  of  the  seismic  load  need  not  neces¬ 
sarily  be  capacitive,  for  it  can  as  well  be  inductive. 


2.  PERFORMANCE  EQUATiONS 

From  the  equivalent  circuit  presented  in  Figure  3,  it  is  possible  to  develop  all  the  relationships 
required  to  characterize  the  performance  of  the  transducer  operating  into  a  seismic  load.  The  quanti¬ 
ties  of  interest  are  the  force  generated  against  the  seismic  medium,  Fm  .  power  radiated  into  the 
medium:  power  conversion  efficiency,  and  information  bandwidth.  This  effort  is  not  concerned  with 


Elactrtal 


Medianical 


1 


1 

Crystal  Stack  Machanical  Resonator  Seismic  Load 


FIGURE  3.  EQUIVALENT  CIRCUIT  REPRESENTATION,  RADC  TRANSDUCER 

DRIVING  SEISMIC  LOAD 

the  design  and  performance  of  the  drive  amplifier  and  associated  elements  since  this  is  considered 
straightforward.  Such  equipment  assembled  to  perform  the  various  tests  will  be  described  later. 


The  transfer  function  relating  the  output  force,  Fm  ,  to  the  input  drive  voltage,  V,  is  given  in  the 
following  expression  where  s  is  the  Laplacian  opera' or  and  the  variables  are  understood  to  be  s- 
domain  variables: 


where: 


hC 


(—£l 


°\Cc+C 


+  C>) 


LCcs2(RCs  +  1) 


V 


/is*  +  Bi 3  4-  Cs 2  +  Ds  +  1 


LCcCCrCpRRr 
Cc  +  C  +  Cr 


(3) 


B=L 


:cCc,r  (i  4-) 


+  CcCd(CR  CrAr)  +  CCrCp(R  +  Rr) 


Cc  +  C  +  Cr 


C  =  L 


CcCR 


Rp(Cc+C  +  Cr ) 


+ 


Cc  TC  +  Cf 


CcCCrRRr 
+  — - — r- 
Cr  +  C  +  Cr 
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Z_  CcCR  +  Cr{CcRr  +  CR  +  CR, ) 

~RP+  Cc  +  C  +  Cr 

This  is  a  rather  formidable  expression,  but,  fortunately,  because  of  the  range  of  numerical 
values  employed  and  in  the  vicinity  of  resonance,  the  transfer  function  is  given  to  a  close  approxi¬ 
mation  by 

hC  ( - 9: - }LCcs2 

Fm_  °\Cc  +  C  +  CrJ  , 


hC0 1 - — - ]LCcs2 

Fm  „  _ °\Cc  +  C  +  CrJ  _  (4 

F  \ro(Q±o  1 ,  +  \L+CCr(R±M]s+l 
Cc  +C  +  Cr  J  [_/?p  Cc  +  C  +  cr_ 

Moreover,  since  the  transducer  must  operate  at  resonance  to  couple  efficiently  to  the  seismic 
medium,  the  magnitude  of  Fm  can  be  determined  from  Eq.  (4)  by  substituting  everywhere  /«  for  s 
and  setting  u>  at  the  value  required  for  resonance.  This  becomes,  for  a  high  Q  resonance, 


'°v{cc+c  +  cr) 


LCcju) 


L_  CCr(R  -1  Rr) 
Rp  Cc  +C  +  Cr 


where  co,  the  resonant  radian  frequency,  is  given  by 


f  Cr(Cc+0 
'  Ce+C  +  C,  p_ 


Note  that  Fm  and  V  are  no  longer  ^-domain  variables  in  Eq.  (5).  The  power  radiated  into  the  seis¬ 
mic  load  is  given  as 

(?) 

Rp 

where  |Fm|  is  given  in  rms  units. 

To  obtain  power  conversion  efficiency,  it  is  necessary  to  determine  the  losses  in  the  transducer 
that  result  from  the  two  lossy  elements  R  and  Rr.  These  two  losses  are  given  by 


Floss 


-fm 


Lu>/  Rr 


The  relationship  between  Fc  and  Fm  is  rather  complex,  but,  for  typical  numerical  values  employed 
here,  sufficient  accuracy  is  achieved  by  setting  |FC|  =  (0.84)|Fm{ .  The  above  equation  can  then  be 


written 


F\oss  ~ 


ifj2  ' 


~(0.706)(/?Cto) 

R 


Rr\L<oJ  J 


As  before,  |Fmj  is  given  in  rms  units.  The  efficiency  of  power  conversion,  ep ,  is  defined  as  follows: 
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£p  = 


rad 


■Prad  "F  -^loss 


(9) 


The  information  bandwidth  is  defined  as  the  bandwidth  at  the  half  power  points  (— 3-dB  points). 
Tliis  is  given  as 


BW  = 


L_  CCr(R  +  Rr ) 


L RP 


Cc+C  +  Cr 


00 
27 r 


,  Hz 


(10) 


As  verified  by  laboratory  tests,  numerical  values  for  all  quantities  are  presented  in  a  later  sec¬ 
tion  of  this  report,  and  performance  calculations  are  presented  in  the  form  of  curves  for  arbitrary 
seismic  load. 

3.  IMPEDANCE  MATCHING 

For  the  applications  envisioned  here,  it  is  desirable  to  minimize  power  supply  requirements 
and  the  weight  and  bulk  volume  of  the  equipment.  Thus  it  becomes  important  to  examine  the 
impedance  matching  problem  since  it  i.,  well  known  that  maximum  power  transfer  occurs  when  the 
source  and  load  impedances  are  matched.  If  we  take  the  deriwrive  of  PJ3a  as  g'.;en  in  Eq.  (1)  with 
respect  to  P.D  and  set  'he  result  to  zero,  we  obtain  the  expression  for  the  value  of  Rp  required  to 
maximize  Prad.  This  is 


R 


L(Cc  +C  +  Cr) 
Popt  CCr(R+Rr ) 


(ID 


Although  not  presented  here,  it  can  be  shown  that  this  value  of  Rp  also  results  in  the  correct  conju¬ 
gate  complex  impedance  match  between  the  seismic  load  and  the  source  transducer,  the  only  other 
requirement  being  that  the  resonant  radian  frequency,  co,  be  that  as  given  by  Eq  (6).  From  this 
fact  it  is  cicar  that  the  resonant  frequency  is  sensitive  to  the  reactive  part  of  the  seismic  load.  As 
shown  by  Miller  and  Pursey(9)  and  as  verified  by  work  under  this  contract,  the  reactive  part  of  the 
radiation  impedance  is  typically  capacitive  where  the  dimensions  of  the  source  are  small  compared 
to  the  shear  wavelength  in  the  medium.  If,  under  other  circumstances,  the  reactive  part  happened 
to  be  inductive,  the  equations  could  readily  be  rewritten.  Then  the  attachment  of  the  transducer  to 
the  seismic  load  would  result  in  an  increase  in  the  resonant  frequency  as  compared  to  the  isolated 
transducer  resonance,  rather  than  in  the  decrease  experienced  with  capacitive-type  loads. 

Miller  and  Pursey(9)  also  show  that  the  magnitude  of  the  radiation  impedance  becomes  greater 
as  the  diameter  of  the  interface  drive  disk  is  increased.  This  implies  that  use  might  be  made  of  area 
transforming  “horns”  between  the  transducer  and  the  seismic  medium  as  a  means  of  obtaining  bet¬ 
ter  impedance  matching.  From  a  ;  'actical  standpoint  in  the  case  at  hand,  this  could  be  accomplished 
only  in  the  direction  of  increasing  the  area  at  the  seismic  medium  end  of  the  horn,  and  this  would 
rapidly  lead  to  very  heavy  and  buiky  horns  with  considerably  aggravated  mounting  problems. 

4.  SIMPLIFIED  EQUIVALENT  CIRCUIT  AND  PERFORMANCE  EQUATIONS 

For  the  range  of  numerical  values  encountered  here  and  in  the  vicinity  of  resonance,  the  equi¬ 
valent  circuit  given  in  Figure  3  can  be  replaced  with  the  simplified  equivalent  circuit  presented  in 
Figure  4. 


f| 


if 

rfi 


if 


1 


w.«£L±5d^.H, 


From  a  conceptual  point  of  view,  this  simplified  model  is  more  satisfying,  and  performance  calcula¬ 
tions  are  less  time  consuming.  The  overall  results  are  identical  with  those  obtained  with  the  more 
complete  expressions.  It  should  be  noted  that  the  facto"  K  is  derived  from  the  comprehensive  cir¬ 
cuit  and  is  given  as 

C, 

*  =  (19) 

5.  THE  TRANSDUCER  AS  A  RECEIVER 

In  all  that  has  gone  before,  the  RADC  transducer  has  been  viewed  only  as  a  source  or  transmit¬ 
ter  for  the  purpose  of  radiating  energy  into  the  seismic  medium.  To  be  sure,  this  is  the  matter  of 
most  direct  concern  here,  since  a  variety  of  small,  low-power,  highly  sensitive  instruments  is  available 
for  the  detection  of  seismic  vibrations.  Nevertheless,  the  RADC  transducer  can  be  employed  as  a 
receiving  device  of  great  sensitivity,  and  it  is  certainly  conceivable  and  manifestly  convenient  that  a 
single  such  unit  at  a  given  site  be  used  for  both  transmitting  and  receiving.  The  equivalent  circuit  for 
the  unit  operating  as  a  receiver  is  as  shown  here: 


Cc\fKhC0: 1 


RP  h  _1F 


The  seismic  disturbance  is  viewed  here  as  a  force  input,  Fs,  while  the  output  is  a  voltage,  Va ,  a 
high  impedance  de.ector  implied.  All  other  circuit  elements  are  as  shown  in  Figure  4.  The  5 -domain 
transfer  function  in  the  vicinity  of  resonance  is,  approximately, 


/ LKhC'\  „  „ 

\  Rpn  js&sCpS+V _ 

KG, +  C-+CV  ♦*(£+£)< 


where- 


H+fPfFCoCc 
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Or,  at  resonance,  substituting  ;cu  everywhere  for  s,  there  results 


V0  ‘ KhC'\(  R*  \  ,  „  _ 
F, ,*  V  n  )(*,  (1  +  '“W 


In  this  form,  this  is  quite  similar  to  Eq.  ( 1 2),  and  it  is  s-oted  that  high  sensitivity  is  predicated  upon 
high  Q  resonance.  In  this  case,  the  resonant  frequency  contains  an  effect  due  to  the  capacitor  across 
the  crystal  terminals,  nC0.  However,  this  can  be  tuned  out  with  an  inductor,  and  the  sensitivity  can 
be  increased  several  decibels  depending  on  the  Q  of  the  inductor.  It  should  be  clear  that  the  resonant 
frequency  of  the  receiver  unit  be  within  a  iew  hertz  of  the  resonant  frequency  of  the  transmitter 
unit  in  order  to  establish  effective  communications. 

6.  LABORATORY  TESTS 

a.  Piezoelectric  Crystal 

Laboratory  tests  on  the  transducer  were  devised  for  the  purpose  of  verifying  the  form  of 
the  mathematical  model  (or  the  equivalent  circuit)  and  quantifying  numerical  values  of  the  various 
elements.  The  first  set  of  tests  was  designed  to  determine  the  values  of  C0 ,  h,  and  Cc  for  the  crystal 
stack.  These  quantities  are  defined  in  the  Appendix.  The  crystal  stack  was  driven  over  the  frequency 
range  of  300  to  400  Hz  with  an  air  load  on  one  face.  This,  in  effect,  short  circuits  the  mechanical 
resonator  and  seismic  load  in  Figure  3.  This  forces  Fc  =  0  in  Eq.  (28),  Appendix,  and,  together  with 
Eq.  (27)  and  considerable  manipulation,  it  is  found  that 


{<  )t+? 


where 


C'  =  — 
juv 

Thus,  the  test  requires  the  measurement  of  V,  u,  i,  and  r  (the  velocity  on  the  air-loaded  face).  Also 
needed  is  the  value  of  n  =  €  and  the  calculated  value  of  k  =  4.58  X  10!  0  N-m  for  the  PZT-4  crystal 
stack  spring  rate.  From  these  data,  values  for  C0  and  h  were  then  calculated.  On  the  following 
page,  the  results  for  each  of  the  transducers  are  tabulated  and  compared  to  overall  theoretical  values 
calculated  from  manufacturers’  data.  It  is  felt  that  these  results,  which  required  veiy  careful  mea¬ 
surements,  demonstrate  the  adequacy  of  the  assumed  model  and  provide  good  numerical  values  for 
purposes  of  calculations. 

b.  Mechanical  Resonator 


The  next  test  was  designed  to  verify  and  quantify  the  mechanical  resonator  portion  of  the 
model.  At  first  glance,  this  appeared  elementary;  however,  the  spring  rate  of  the  main  resonator 
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Quantity 


Trans.  No.  30  I  Trans.  No.  31  I  Theoretical 


Co.lt  F 

0.00227 

0.00213 

0.00204 

h,  V/m  X  109 

2.70 

2.94 

3.15 

hC0,  N/V 

6.12 

6.26 

6.42 

nC0,liF 

0.0136 

0.0128 

0.0122 

Cc,  m/N  X  10‘10 

2.05 

2.27 

2.35 

••  p 

0.0212 

0.0215 

0.0219 

spring  was  found  to  be  of  the  order  of  106  lb/in.  such  that  there  was  simply  notiling  “stiff’  enough 
against  which  this  spring-mass  system  might  be  referenced  for  measurement  purposes  without  also 
affecting  the  resonance  and  Q  of  the  system.  The  solution  was  to  mount  the  two  transducers  back- 
to-back  and  drive  them  in  parallel.  In  this  arrangement,  assuming  that  the  two  units  are  well  matched 
in  dynamic  properties,  the  interface  between  the  t?»o  units  appears  as  a  nearly  infinitely  “stiff’  refer¬ 
ence,  and  the  resulting  measurements  of  resonant  frequency  and  Q  are  representative  of  the  isolated 
mechanical  resonator.  Accordingly,  this  test  was  performed  with  typical  results  presented  in 
Figure  5.  This  is  a  frequency  response  curve  presenting  the  ratio  of  output  force  to  input  drive 


Frsquency,  •'-Hz 


FIGURE  5.  FREQUENCY  RESPONSE  RADC  TRANSDUCER  INFINITELY 
STIFF  LOAD  (BACK-TO-BACK) 
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voltage  over  the  frequency  range  of  250  to  800  Hz.  Output  force  was  determined  as  the  product  of 
mechanical  resonator  mass  and  the  measured  acceleration  of  the  mass.  Input  drive  voltage  and  fre¬ 
quency  were  measured,  and,  in  addition,  an  accelerometer  was  mounted  on  the  interface  plate 
between  the  two  transducers  to  moni'or  the  quality  of  the  dynamic  match  between  the  two  units. 
This  signal  was  found  to  be  more  than  40  dB  less  than  the  mass  accelerations,  and  so  the  match 
appeared  adequate  for  purposes  here.  T!.is  was  further  verified  by  the  fact  that  the  measured  accel¬ 
erations  on  each  of  the  two  masses  were  very  nearly  equal.  With  calculated  values  of  L  and  Cr  and 
with  already  measured  values  of  hC0  and  Cc,  the  results  obtained  from  Figure  5  were  sufficient  to 
determine  values  for  C  and  ( R  +  Rr)  [see  Fqs.  (5)  and  (6)  with  Rp=°°  and  Cp  =  0] .  R  and  Rr  were 
then  apportioned  to  the  two  springs  in  inverse  proportion  to  the  respective  spring  rates. 

c.  Tabulated  Numerical  Values 


From  the  preceding  laboratory  tests  and  by  calculation  of  the  more  certain  quantities,  the 
following  numerical  values  were  assigned  for  the  elements  shown  in  the  equivalent  circuit  model  of 
Figure  3.  The  values  have  been  averaged  for  each  of  the  two  RADC  transducers  and  are  given  in 
both  MKS  and  convenient  engineering  units. 

nC0  -  0.0132  uF 
hC0  =  6.2  A  sec/m 
hCQ  =  6.2  N/V 
Cc  =  2.1 2  X  l«r10  m/N 
Cr  =  233  X  10*10  m/N 
C  =  42.8  a  10'10  m/N 
R  =  690  N  sec/m 
Rr  -  126  N  sec/m 
L  =  56.5  N  sec?/m 


=  0.158  A  sec/in. 

=  1.4  lb /V 

=  0.037  X  10'6  in./lb 
=  4.07  X  10'6  in./lb 
=  0.748  X  10‘6  in./lb 
=  3.94  lb  sec/in. 

=  0.72  lb  sec/in. 

=  0.323  lb  sec2  /in. 


If  one  chooses  to  employ  the  simplified  equivalent  circuit  model  of  Figure  4  and  the  related  equa¬ 
tions,  the  following  values  are  used: 


nCQ  =0.0132  mF 
hC0  =  6.2  A  sec/m 
Cc  =  2.12  X  10'10  m/N 
L  -  56.5  N  sec2/m 
C*  =  35.5  X  iO'10  m/N 
R*=  19.4  X  106  Nsec/m 
K  =  0.837 


=  0.158  A  sec/in. 

=  0.037  X  1CT6  in./lb 
=  0.323  lb  sec2 /in. 

=  0.62  X  10’6  in./lb 
=  111,000  lb  sec/in. 


d.  Performance  Calculations 


Using  the  numerical  values  cited  above,  performance  calculations  were  made,  and  ihe 
results  are  shown  plotted  in  Figures  6  and  7.  The  first  curve  presents  the  radiated  power  output  for 
various  drive  voltages  as  a  function  of  the  equivalent  parallel  seismic  load  resistance,  Rp ,  while  the 
second  curve  indicates  the  variation  of  power  conversion  efficiency  and  information  bandwidth  with 
this  same  parameter.  The  curves  were  calculated  for  an  assumed  nominal  resonant  frequency  of 
334  Hz.  Several  interesting  points  are  seen  at  once  by  examination  of  these  curves.  First,  the  power 
radiating  capability  of  the  transducer  is  strongly  load  dependent  and  increases  with  increasing  Rp  up 
to  the  matched  load  point.  Higher  values  of  Rp  are  associated  with  hard,  competent  rock. 
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Consequently,  the  unit  w,m  clearly  designed  for  operation  Into  such  media,  Maximum  power  opera* 
tion  at  the  matched  loud  point  cannot,  In  fact,  be  realized  because  of  motional  limits  on  the  excur¬ 
sion  of  the  resonant  nu««« ,  such  limits  were  probably  In  turn  derived  from  maximum  operating  stress 
limits.  Rock  media  exhibiting  such  a  high  value  of  /Jy,  "III  ,000  lb  ln,/*cc  are  probably  nonexistent 
in  any  event,  and  so  act ual  maximum  power  operation  In  the  vicinity  of  Rp  ■  40,000  to  60,000  lb  see / 
in,  is  more  realistic,  llccuuse  of  narrow  practical  limits,  there  was  no  attempt  here  to  investigate 
impedance  matching  horn-,  beyond  that  of  providing  the  6-in, -diameter  drive  plate  in  lieu  of  using 
the  3-ln,*diameter  transducer  base.  This  appeared  to  be  a  step  in  the  right  direction  from  the  stand¬ 
point  of  impedance  mi  telling.  Second,  the  transducer  does  indicate  high-power  conversion  efficien¬ 
cies  In  the  range  of  0,8  or  so  for  typical  hard  rock  media  encountered  in  field  experiments  to  be  des¬ 
cribed  later.  And  third,  the  device  exhibits  a  very  limited  Information  bandwidth  of  the  order  of 
10  to  12  Hz,  The  nanow  bandwidth  is  not  surprising  when  the  high  Q  nature  of  the  unit  is  con¬ 
sidered,  Moieovcr,  the  high  Q  characteristic  is  not  incidental;  it  has,  in  fact,  been  employed  to  mul¬ 
tiply  the  reaction  of  the  m  m  element  (125  lb  weight)  acting  against  the  outer  surface  of  the  crystal 
stack  and  thus  facilitating  the  radiation  of  power  into  the  seismic  medium,  To  make  the  unit  broad- 
banded  with  the  same  power  radiating  capability  would  require  the  use  of  a  reaction  mass  (non- 
resonant)  weighing  10,000  lb  or  more.  Not  only  is  this  an  absurdity,  but  such  a  mass  of  steel  or  lead 
would  doubtless  corrupi  the  frequency  range  of  interest  with  a  whole  spectrum  of  internal  propa- 
gational  modes;  i,c,,  it  would  not  appear  as  a  lumped  reactive  mass  element  at  the  crystal  stack  face. 
Another  solution  would  be  to  load  both  faces  of  the  ciystal  stack  against  the  rock  medium  by 
mounting  the  unit  in  a  crevice  between  two  rock  faces,  for  example;  again,  this  would  not  be  too 
practical  an  approach  for  general  use. 


o.  Powor  Transfer  (Motor-Generator  Mode) 


Since  little  power  was  consumed  in  the  back-to-back  tests  reported  previously,  it  was 
decided  to  attempt  to  transfer  power  in  a  motor-generator  mode  in  which  the  two  transducers  were 
mounted  back-to-back  with  the  drive  applied  to  one  transducer  and  the  other  loaded  with  an  electri¬ 
cal  resistor  across  the  crystal  terminals.  It  was  thought  that,  in  this  way,  maximum  power  radiating 
capabilities  could  be  determined  as  well  as  efficiency  of  power  conversion.  However,  these  tests  did 
not  prove  particularly  satisfying  for  two  reasons,  First,  extending  the  simplified  equivalent  circuit 
of  Figure  4  to  represent  the  motor-generator  mode,  there  results  the  following  equivalent  circuit; 
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The  input  power  is  given  by  Vi  cos  6  where  6  is  the  phase  angle  while  the  output  power  is  J V^Ri 
where  all  quantities  are  given  in  rms  units.  The  problem  arises  in  that  this  circuit  does  not  behave 
simply  as  a  single  transducer  driving  into  a  seismic  load  as  might  be  suspected  by  inspection.  For 
one  thing,  the  Q  at  resonance  is  lower,  and  it  might  be  expected  that  internal  losses  would  be  higher 
by  a  factor  of  two.  Also  the  electrical  load,  Ri,  is  not  linearly  related  to  the  seismic  load  element, 

Rp ,  and  in  fact,  it  is  not  possible  to  simulate  a  seismic  load  less  than  65,000  lb  sec/inch.  In  addition, 
it  was  found  that,  because  of  the  nonlinear  behavior  of  the  transducer,  the  current  wave  shape  was 
not  sufficiently  sinusoidal  to  permit  accurate  phase  angle  measurements  for  purposes  of  calculating 
input  power.  Despite  these  conceptual  and  practical  difficulties,  tests  were  run  in  which  output 
power  as  high  as  25  W  was  measured  with  input  power  of  perhaps  50  W  corresponding  to  an  equi¬ 
valent  Rp  of  about  65,000  lb  sec/inch.  By  interpretation,  assuming  the  losses  were  twice  as  high  in 
this  mode,  this  would  correspond  to  an  efficiency  of  about  0.67  for  a  single  transducer  operating 
into  the  same  seismic  load.  In  Figure  7,  it  is  seen  that  the  correct  efficiency  for/?p  =  65,000  lb  sec/ 
in.  is  0.63.  Consequently,  this  exercise  did  serve  to  corroborate  previous  results  in  general,  as  well 
as  to  demonstrate  the  transfer  of  a  reasonable  amount  of  power. 


f.  Transducer  Nonlinearity 


The  problem  of  nonlinear  behavior  of  the  transducer  was  just  mentioned;  thus,  it  seems  in 
order  to  further  discuss  this  matter  since  the  ramifications  of  this  difficulty  were  felt  throughout  the 
remainder  of  this  work.  The  initial  symptoms  of  this  nonlinearity  were  (1 )  a  definite  downshift  in 
the  transducer  resonant  frequency  with  increasing  output  level,  making  it  extremely  difficult  because 
of  the  high  Q  to  tune  onto  resonance  with  an  oscillator  and  (2)  the  observed  nonsinusoidal  depar¬ 
tures  of  the  drive  current  wave  shape.  Later,  when  the  drive  plate  was  instrumented  with  acceler¬ 
ometers,  the  nonsinusoidal  wave  shape  was  noted  again.  In  all  cases,  the  nonlinearity  increased  with 
output  level.  Interestingly,  the  measured  acceleration  of  the  resonant  mass  did  not  deviate  as 
markedly  from  sinusoidal  shape,  apparently  because  of  the  filtering  effect  of  the  high  Q  resonator. 

In  diagnosing  this  problem,  it  was  found  that  the  nonlinearity  occurred  only  on  the  half  cycle  when 
the  operating  forces  were  in  a  direction  to  relieve  the  load  applied  by  the  preload  rod.  The  main 
spring  in  effect  becomes  softer  as  the  total  load  is  decreased.  Tightening  of  the  preload  rod  improved 
the  nonlinearity,  but  this  effect  could  not  be  pursued  because  of  the  already  high  working  stresses  in 
the  rod  (SAE  4140,  heat  treated).  As  noted  earlier,  the  crystal  stack,  main  spring,  and  mass  are  all 
coupled  solely  through  this  preload,  and  it  is  believed  that  this  is  a  design  deficiency  from  which  the 
nonlinearity  originates.  This  could  be  remedied  through  proper  redesign.  At  any  rate,  the  problem 
did  affect  all  succeeding  tests  where  drive  plate  acceleration  data  were  measured  for  the  purpose  of 
calculating  radiated  power  and  driving  point  impedance,  since  it  was  not  possible  to  measure  phase 
angle  with  any  reasonable  degree  of  precision. 

g.  Driving  Concrete  Pier  Load 

1  o  provide  some  foreknowledge  of  transducer  performance  operating  into  a  seismic  load 
and  to  provide,  simultaneously,  a  test-bed  for  finalizing  the  experimental  equipment  design,  the 
transducer  was  attached  to  the  concrete  floor  in  the  laboratory',  directly  above  a  heavy  concrete 
floor  beam  and  an  18-in.-diameter  pier  extending  some  36  ft  down  to  limestone  layers  in  the  earth. 

A  heavy  concrete  anchor  bolt  imbedded  in  the  floor  was  drilled  and  tapped  to  accept  the  preload 
rod  for  purposes  of  positive  attachment.  A  5-in.-diamcter  drive  plate  was  employed  in  this  instance, 
principally  to  provide  a  pad  for  mounting  an  accelerometer.  Thus,  the  transducer  was  affixed  to  the 
concrete  at  a  single-point  tic  using  the  preload  rod.  From  this  experience,  during  which  the  crystal 
stack  on  one  of  the  transducers  was  fractured  because  of  inadvertant  loss  of  crystal  preload  during 


the  mounting  procedure,  it  was  decided  to  employ  a  different  mounting  technique  for  the  field 
experiments.  A  6-in.-diameter  drive  plate  was  designed,  providing  a  means  for  anchoring  the  preload 
rod  independently  of  attachment  to  the  rock  medium.  This  in  turn  required  the  use  of  three  smaller 
anchor  bolts  to  attach  the  transducer  to  the  rock.  By  this  method,  the  transducer  assembly  was 
always  under  preload  and  not  subject  to  accidental  damage  during  handling,  mounting,  etc. 

Both  RADC  transducers  were  operated  into  the  concrete  pier  load.  A  typical  frequency 
response  wave  for  one  of  the  units  is  presented  in  Figure  8.  This  gives  the  amplitude  ratio  of  out¬ 
put  force  to  drive  voltage,  and  the  loss  in  Q  due  to  loading  is  seen  by  comparing  this  to  the  infinite 
load  frequency  response  of  Figure  5.  Note  also  the  lower  resonant  frequency.  The  driving  point 
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FIGURE  8.  FREQUENCY  RESPONSE  RADC  TRANSDUCER  CONCRETE  PIER  LOAD 

impedance  at  th~  *Tansducer  concrete  floor  interface  was  calculated  from  measured  values  of  inter¬ 
face  force,  Fm  ,  and  interface  acceleration  (converted  to  velocity).  This  required  phase  angle  mea¬ 
surement,  and  precision  was  therefore  limited  by  nonlinear  behavior  as  discussed  above,  although 
these  measurements  were  made  at  low  drive  levels  with  correspondingly  reduced  nonlinear  effects. 
The  average  value  of  driving  point  impedance,  Zg.  determined  from  a  number  of  individual  mea¬ 
surements  was,  at  a  resonant  frequency  of  337  Hz : 
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Zg  =  4000  — /II  ,500  lb  sec/in. 

This  corresponds  to 

Rp  =  37,000  ib  sec/in. 

Cp  =0.037  X  1  O'6  in./lb 

Moreover,  it  was  found  that  the  paralled  representation  with  the  values  given  here  was  quite  good 
for  this  concrete  pier  load  over  the  range  of  200  to  800  Hz.  Plugging  these  values  of  Rp  and  Cp 
back  into  Eqs.  (5)  and  (6),  or  alternatively  Eqs.  (12)  and  (13),  yields  values  of  Fm/V  and  a?  at 
resonance  very  close  to  those  given  in  the  experimental  data.  This  further  confirms  the  form  of  the 
mathematical  model  and  the  numerical  values  previously  determined. 

Power  radiated  into  the  load  is  given  by  Fm i>  cos  $  where  v  is  the  velocity  of  the  drive 
plate  and  0  is  the  phase  angle  between  Fm  arid  >\  Fm  and  r  are  given  in  rms  units.  Radiated  power 
was  calculated  from  the  experimerital  data  using  this  expression  and  compared  to  the  radiated  power 
given  in  Figure  6  using  Rp  =  37,000  lb  sec/in.  and  the  drive  voltage  level.  These  comparisons  were 
not  consistent  over  a  number  of  separate  measurements  evidently  because  of  the  lack  of  precision  in 
phase  angle  measurement.  It  was  concluded  that  the  best  estimate  would  result  by  using  Figure  6 
with  the  measured  value  of  .Rp.  From  this  it  is  seen  that  power  c  .'about  43  W  was  radiated  into  the 
concrete  pier  for  a  drive  voltage  of  1500  V  rms.  This  same  technique  was  later  used  in  the  field 
experiments  to  estimate  radiated  power.  i 


The  electrical  input  impedance,  Z;,  at  the  crystal  stack  terminals  was  found,  for  the  same 
condition  of  resonance  as  above,  to  be 

Zi  =  9750-/  16,900  S2 

The  corresponding  input  power  for  V  =  1 500  V  rms  is  57.5  W,  which  yields  a  power  conversion 
efficiency  of  0.75,  in  very  ciose  agreement  with  the  result  given  in  Figure  6.  Note  also. that  the 
drive  force,  Fm ,  was  3810  lb  rms  (5390  ib  peak). 

* 

These  concrete  pier  tests  also  facilitated  the  completion  of  the  experimental  equipment 
buildup,  the  final  form  of  which  is  described  in  the  next  section.  : 


19 


SECTION  ill 


EXPERIMENTAL  EQUIPMENT 


1.  GENERAL  ARRANGEMENT 

Development  of  instrumentation  and  equipment  was  not  one  of  the  contract  objectives.  How¬ 
ever,  a  buildup  of  equipment  was  required  to  perform  the  various  laboratory  and  field  tests,  and  it 
might  be  instructive  to  review  briefly  those  arrangements.  For  the  field  experiments,  two  general 
types  of  tests  were  planned: 

(1 )  Local  transmitter  tests  employing  CW  excitation  and  instrumentation  to  measure  resonant 
mass  and  drive  plate  accelerations,  from  which,  together  with  other  measured  data,  the 
driving  point  impedance,  Zg,  could  be  calculated.  This  in  turn  permitted  an  estimate  of 
radiated  power,  conversion  efficiency,  force  output,  and  bandwidth  at  the  transmitter 
installation. 

(2)  Seismic  propagation  and  communications  experiments,  employing  one  RADC  transducer 
as  the  transmitter  and  the  other  as  a  receiver,  with  operation  in  a  gated  pulse  mode. 

Because  of  the  narrow  information  bandwidth  and  poor  control  of  the  transducer-installed 
resonant  frequency,  it  appeared  that  little  would  be  gained  by  consideration  of  modula¬ 
tion  methods  other  than  simple  pulse  width  or  pulse  code.  Thus,  it  was  planned  to  trans¬ 
mit  in  the  gated  pulse  mode  with  capability  of  varying  the  on-time  of  the  pulse.  Means  for 
remote  triggering  of  the  transmitter  from  the  receiver  site,  via  wire  line,  was  provided  so 
as  to  facilitate  propagation  velocity  measurements.  Provision  was  also  made  to  employ  a 
sensitive  accelerometer  as  a  seismic  receiving  unit. 

Block  diagrams  showing  the  equipment  arrangements  for  the  local  transmitter  tests  and  for  the 
transmitter  and  receiver  unit  propagation  tests  are  presented  in  Figures  9,  10,  and  1 1 ,  respectively. 
Corresponding  photographs  of  the  equipment  packages  are  presented  in  Figures  12,  13,  and  14, 
respectively.  Standard  laboratory  instrumentation  was  employed  wherever  possible  though  this  did 
result  in  sporadic  operational  difficulties  in  the  iron  mine  tests  due  to  the  lugh  humidity  environ¬ 
ment.  Specia1  equipment  design  and  fabrication  were  required  in  only  two  areas:  (1)  the  high-volt¬ 
age  drive  transformer  and  (2)  the  unit  identified  as  the  control  circuit.  A  description  of  these  two 
items  follows. 

2.  DRIVE  TRANSFORMER 

This  unit  was  required  to  match  the  dual,  ciass  A.  tube-type  audio  power  amplifier  to  the  trans¬ 
ducer  crystal  stack.  The  dual  amplifiers  were  connected  to  provide  an  output  impedance  of  32  £2  and 
a  peak  voltage  of  about  57  V  rms.  while  delivering  a  maximum  power  of  100  W.  Input  power 
required  from  the  1 10-VAC  line  was  200  W.  The  crystal  stack,  on  the  other  hand,  required  that 
1500  V  rms  be  delivered  into  an  impedance  of  the  order  of  (9750  — /  16,900  £2)  at  the  nominal  reso¬ 
nant  frequency  of  about  335  Hz.  A  special  transformer  was  ordered  with  a  27:1  step-up  ratio,  and 
a  10-henry  choke  was  connected  across  the  secondary  to  correct  the  pool  power  factor  of  the  load. 
This  resulted  in  a  fairly  good  match  with  the  load  appearing  as  a  nearly  resistive  lead  of  about 
39,000  a. 
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Note:  All  equipment  battery  powered. 


3.  CONTROL  CIRCUIT 


In  the  laboratory  tests,  there  was  considerable  difficulty  in  tuning  the  transducer  to  resonance. 
Further,  since  the  resonance  is  shifted  by  the  seismic  load  reactance  at  each  installation,  there  could 
be  no  predetermined  drive  frequency.  Accordingly,  it  was  decided  to  close  the  loop  around  th ; 
transducer  and  allow  it  to  oscillate  at  its  own  self-resonance  and  thus  generate  maximum  power 
radiating  capability.  This  was  accomplished  by  installing  a  velocity  sensor  on  the  RADC  transducer 
resonant  mass  and  feeding  this  signal  back  to  the  drive  amplifier  in  positive  feedback  manner.  The 
control  circuit,  shown  schematically  in  Figure  15,  was  devised  to  accomplish  this  purpose.  Since  it 


FIGURE  15.  CONTROL  BOX-SIMPLIFIED  SCHEMATIC 

was  desired  to  operate  the  source  transducer  in  a  gated  lulse  mode,  provisions  were  made  to  close 
and  open  the  feedback  loop  with  a  solid-state  switch  controlled  by  a  one-shot.  The  one-shot  in 
turn  was  fired  by  the  remote  trigger  wire  line.  Further  refinements  were  incorporated  in  the  circuit 
to  permit  very  fast  tum-on  of  the  transmitter  limited  only  by  the  inherent  Q  of  the  mechanical 
resonance.  One  of  these  consisted  of  a  variable  gain  element,  a  multiplier,  which  was  configured  to 
provide  very  high  loop  gain  at  the  initiation  of  a  pulse  but  with  gain  rapidly  decreased  by  the  build¬ 
up  of  the  average  loop  drive  level.  Thus,  the  final  drive  level  was  controlled  but  without  inhibiting 
the  initial  buildup  time.  Also  a  start  “kick”  pulse  was  provided  to  bring  the  velocity  sensor  signal 
up  to  the  threshold  of  response.  The  phase  rhift  adjustment  was  found  useful  in  permitting  reso¬ 
nance  to  occur  directly  at  the  point  of  maximum  amplitude  excursion.  The  control  circuit  was 
successfully  employed  in  the  concrete  pier  experiments  and  later  in  the  field;  however,  each  separ¬ 
ate  installation  required  appropriate  tune-up  of  gains,  phase  shift,  etc. 


SECTION  IV 


LIMESTONE  QUARRY  EXPERIMENTS 


1.  TEST  SITE  DESCRIPTION 

Seismic  communications  tests  were  initiated  in  an  open-formation  limestone  quarry  about 
35  mi  northeast  of  San  Antonio  near  New  Braunfels,  Texas.  The  quarry  is  owned  and  operated  by 
U.S.  Gypsum  Co.  These  tests  were  performed  over  a  period  covering  the  first  2  weeks  in  April  1972. 
A  plan  view  of  the  test  site  and  a  typical  cross  section  are  presented  in  Figure  1 6.  The  limestone  is 
identified  as  Edwards  Limestone  and  is  considered  a  relatively  hard,  competent  rock  generally, 
though  varying  in  places  to  relatively  soft  and  spongy  in  character.  As  seen  in  the  figure,  one  receiver 
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site  was  employed,  while  the  transmitter  was  installed  at  two  different  sites.  The  RADC  transducer 
installation  at  the  receiver  site  is  shown  in  Figure  17,  while  Figure  18  is  a  photograph  of  a  typical 
mounting  “pad”.  The  anchor  bolts  are  of  the  split-ring  type  and  are  3/4  in.  in  diameter  and  6  in. 


long.  Quick-setting  grout  was  applied  under  the  drive  plate  to  provide  uniform  coupling  to  the  rock 
across  the  full  face  of  the  plate.  Figure  19  is  a  photograph  of  the  transducer  installation  at  trans¬ 
mitter  site  No.  1 ,  showing  plainly  the  layered  limestone  structure.  The  other  two  transducer  instal¬ 
lations  appeared  very  similar  to  this,  and,  in  each  case,  the  transducer  was  affixed  to  what  was 
believed  to  be  the  identical  limestone  layer.  At  each  site,  an  attempt  was  made  to  select  an  attach¬ 
ment  point  on  the  rock  that  appeared  to  be  free  from  surface  cracks  and  that  evidenced  hard,  com¬ 
petent  character  in  response  to  hammer  blows. 


a.  Propagation  Path:  450  Ft 

So-called  “local  transmitter”  tests  were  performed  at  both  the  initial  transducer  installa¬ 
tions,  identified  later  as  the  receiver  site,  and  at  transmitter  site  No.  1 .  It  was  found  that  little 
power  could  be  radiated  at  the  receiver  site.  This  was  corroborated  by  driving  point  impedance 
measurements  yielding 


Rp  =  8370  lb  sec/in. 

Cp  =0.165  X  10'6  in./lb 

at  a  resonant  frequency  of  317  Hz.  Measured  force,  Fm ,  was  637  lb  rms  at  maximum  drive  voltage 
capability  of  950  V  rms.  The  radiated  power  worked  out  to  be  about  5  W.  At  transmitter  site 
No.  I ,  the  results  v  ere  much  superior,  with  driving  point  impedance  measurements  giving 

Rp  =  17,800  lb  sec/in. 

Cp  =0.143  X  10'6  in./lb 

at  a  resonant  frequency  of  315  Hz.  Measured  force  was  1975  lb  rms  at  1500  V  rms,  and  radiated 
power  was  about  25  W.  No  explanation  for  this  discrepancy  in  driving  point  impedance  at  different 
points  in  the  same  rock  layer  was  found.  However,  it  is  surmised  that  subsurface  cracks  in  the  rock 
formation,  perhaps  even  several  feet  below  the  exposed  surface,  could  well  result  in  variations  of 
this  magnitude.  At  any  rate,  because  of  these  findings,  signal  transmissions  were  made  only  from 
the  latter  site  and  received  ♦he  former,  thus  the  identifications  as  given  above. 

The  transmitter  unit  was  triggered  by  wire  line  from  the  receiver  site,  and  a  large  number 
of  pulse  signals  were  recorded  with  signal-to-noise  ratios  greater  than  40  dB  despite  the  presence  of 
heavy  quarry  machinery  operating  about  200  yd  away.  Actually,  acoustic  coupling  due  to  near  air¬ 
craft  flyovers  and  close  human  voices  was  more  of  a  problem.  Figure  20  presents  typical  recorded 
signals.  The  upper  trace  is  the  trigger;  the  second  is  the  acceleration  of  the  transmitter  unit  drive 
plate  (the  seismic  input  signal  in  this  cast);  and  the  third  is  the  signal  received  by  the  receiver  unit. 
The  input  gated  pulse  is  about  0.14  sec  long  and  exhibits  essentially  first-order  exponential  rise  and 
decay  as  determined  by  the  Q  of  system  resonance.  (Note  the  start  “kick”  at  the  trigger  point.)  The 
received  signal  exhibits  second-order-type  onset  just  as  it  might  be  expected  for  a  first-order  system 
to  respond  to  a  first-order  signal  input.  Note  that  the  tail  trails  out  considerably. 

The  velocity  of  propagation  of  this  signal  was  found  by  subtracting  the  time  for  the  input 
signal  to  reach  63  percent  of  final  value  on  the  onset  from  the  time  required  for  the  received  signal 
to  reach  27  percent  of  final  value.  This  is  the  correct  relationship  if  the  two  first-order  time  con¬ 
stants  are  equal.  This  time  difference  was  0.05  sec.  The  propagation  velocity  for  the  450-ft  path 
was  thus  9000  ft/ sec.  Since  this  velocity  appeared  slow  for  the  compressional  wave  velocity  in  hard, 
competent  limestone,  impulse  tests  were  performed  to  shed  some  light  on  propagational  velocities  in 
this  seismic  transmission  path.  Impulses  were  generated  by  sledge  hammer  blows  on  the  rock  face 
within  a  few  feet  of  the  transmitter  unit  installation.  The  drive  plate  accelerometer  signal  was  again 
recorded  as  the  input  signal,  and  received  signals  were  recorded  from  the  RADC  transducer  and  from 
a  sensitive  accelerometer  at  the  receiver  site.  Typical  results  of  these  tests  are  shown  in  Figure  2 1 . 
The  upper  trace  is  the  seismic  input  signal;  the  second  is  the  received  signal  from  the  RADC  trans¬ 
ducer;  and  the  third  is  the  received  signal  from  the  sensitive  accelerometer.  The  RADC  transduce! 
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INPUT  SIGNAL:  TRANSMITTER  SITF  NO.  1 
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RECEIVED  SIGNAL:  RADC  TRANSDUCER 


FIGURE  20.  LIMESTONE  QUARRY-TYPICAL  TRANSMISSION  (450-Ft  Path) 

is  of  course  nam>w  band  with  a  bandwidth  of  perhaps  20  Hz  centered  about  3 1 7  Hz.  The  sensitive 
accelerometer  on  the  hand  is  broadband  over  a  wide  range  leading  up  to  a  resonance  around  800  Hz. 
The  first  arrival  on  both  receiving  units  is  seen  to  come  at  about  0.025  sec  from  the  initial  onset  of 
the  impulse.  This  gives  a  propagation  velocity  of  1 8,000  ft/sec,  which  is  a  reasonable  value  for  the 
compressional  wave  velocity  in  the  medium.  This  confirms  that  compressional  wave  propagation  was 
insignificant  in  “he  RADC  transducer  gated  pulse  tests  and  that  the  dominant  mode  might  have  been 
shear  body  wives,  surface  waves,  or  some  other  form  of  more  complex  waves  as  discussed  by  Owen.n) 
Miller  and  Pur  vey(10)  show  that,  for  a  source  driving  into  an  isotropic  half  space,  the  energy  is  parti¬ 
tioned  approximately  as  follows: 

compressional  wave:  7  percent 
shear  wave:  26  percent 
surface  wave:  66  percent 

This  indicates  that  surface  wave  propagation,  at  least  under  the  ideal  conditions  assumed,  soaks  up 
most  of  the  energy  available.  Accordingly,  it  is  probably  not  too  inaccurate  to  assume  that  the  major 
propagationa!  mode  encountered  here  was  surface  wave,  with  the  wave  simply  propagating  over  the 
open  formation  with  a  wavelength  of  about  30  ft. 

Since  i  e  main  applications  of  interest  here  preclude  surface  wave  propagation,  simply 
from  topological  considerations,  it  was  concluded  that  further  tests  in  an  actual  mine  environment 
would  be  necessary  to  ascertain  the  feasibility  of  the  subject  seismic  communications  method. 
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RECEIVED  SIGNAL:  SENSITIVE  ACCELEROMETER 


FIGURE  21 .  LIMESTONE  QUARRY-IMPULSE  RESPONSE  (450-Fl  Path) 

Accordingly,  arrangements  were  made  to  mount  an  experimental  effort  in  an  operating  mine,  and 
these  efforts  and  the  results  achieved  are  reported  in  the  next  section.  However,  to  gain  seme  fur¬ 
ther  experimental  data  in  the  quarry,  transmission  over  a  longer  distance  was  attempted  with  results 
reported  next. 

b.  Propagation  Path:  900  Ft 

One  RADC  transducer  was  moved  to  transmitter  site  No.  2  as  shown  in  Figure  16,  while 
the  receiver  unit  was  retained  at  the  original  location.  As  before,  “local  transmitter”  tests  were  per¬ 
formed  with  the  driving  point  impedance  measurements  giving 

Rp  =  9100  lb  sec/in. 

Cp  =0.103  X  10'6  in./lb 

at  a  resonant  frequency  of  320  Hz.  Measured  force  was  1020  lb  rms  at  a  maximum  drive  voltage 
capability  of  1375  V  rms.  Radiated  power  was  about  14  W.  Again,  as  at  the  receiver  site,  this 
proved  to  be  a  relatively  “soft”  coupling.  The  loss  in  radiated  power,  the  doubling  of  the  length  of 
the  transmission  path,  and  the  existence  of  what  appeared  to  be  large  fractures  across  the  transmis¬ 
sion  path  all  combined  to  the  extent  that  no  signal  was  received  above  the  noise  level  tor  any  of 
several  types  of  receivers.  Nor  was  it  possible  to  receive  impulse  signals  generated  by  sledge  hammer 
blows.  Thus,  some  idea  of  potential  limitations  to  this  method  of  seismic  communications  was  made 
evident. 
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basement  rock  chemically  identical  to  granite  but  having  finer  grain  structure.  In  consultation  with 
the  mine  people,  locations  were  selected  such  that  no  known  dikes  or  fractures  intercepted  the 
transmission  paths  and  for  which  there  were  no  direct  interconnecting  paths  that  might  support 
significant  surface  wave  propagation.  Other  features  that  had  to  be  considered  included  the  avail¬ 
ability  of  1 10-VAC  power  at  ihe  transmitter  site  and  of  an  air  supply  at  all  sites  to  facilitate  drilling 
of  the  anchor  bolt  holes;  relative  remoteness  from  heavy  machinery,  vent  fans,  etc.;  noninterference 
with  mine  operations;  and  means  for  direct  wire  line  interconnection  between  the  transmitting  and 
receiving  sites.  It  was  not  possible  to  meet  all  of  these  requirements  at  each  site,  nor  was  it  possible 
to  perform  the  full  spectrum  of  tests  that  would  have  been  desirable.  Mine  management  had  author¬ 
ized  experimental  work  during  the  2-week  period  covering  the  second  and  third  weeks  of  May  1972, 
and  it  should  be  noted  that  these  efforts  required  rather  considerable  support  from  mine  personnel 
and  consequent  disruption  of  normal  activities.  The  logistical  and  environmental  difficulties  asso¬ 
ciated  with  mounting  this  sort  of  experimental  work  in  the  far  reaches  of  an  operating  mine,  particu¬ 
larly  involving  the  use  of  delicate  electronic  equipment,  seriously  limited  the  scope  and  extent  of 
work  that  could  be  accomplished  in  a  2-week  period.  However,  the  intent  here  is  not  to  dwell  upon 
how  these  difficulties  limited  the  work  but  rather  to  report  upon  what  was  accomplished,  in  sum, 
the  transmitter  unit  was  installed  on  the  2475  level,  and  signals  were  received,  more  or  less  vertically 
above  that  point,  first  on  the  2275  level  and  then  on  die  1 825  level.  From  directly  observable 
results  in  the  mine,  it  did  not  appear  that  an  attempt  to  transmit  over  a  longer  path  was  justified. 

As  will  be  shown  later,  that  judgment  was  premature  and  incorrect.  At  any  rate,  the  expiration  of 
schedule  time  precluded  further  experimental  work  of  significance.  Summary  descriptions  of  the 
experimental  arrangements  and  results  obtained  follow. 

2.  TEST  RESULTS 

a.  Transmitter  Installation:  2475  Level 

One  RADC  transducer  was  installed  in  the  back  of  d  drift  (roof  of  a  tunnel)  on  the  2475 
level.  A  photograph  of  this  installation  is  presented  in  Figure  23.  Local  transmitter  tests  indicated 
that  good  coupling  to  the  rock  body  was  achieved  with 

Rp  =  25,000  lb  sec/in. 

Cp  =0.15  X  10'6  in./lb 

at  a  resonant  frequency  of  308  Hz.  With  1 500-V  rms  drive,  the  force  was  2570  lb  rms,  and  the  peak 
radiated  power  was  about  30  W.  The  power  amplifier  consumed  about  200  W  from  the  1 10-VAC 
line. 

b.  Receiver  Installation:  2275  Level 

The  other  RADC  transducer  was  installed  in  the  rib  o»  a  drift  (wall  of  a  tunnel)  on  the 
2275  level.  From  mine  survey  records,  the  vertical  distance  to  the  transmitter  was  183  ft  and  the 
horizontal  offset  165  ft.  yielding  a  slant  range  through  the  porphyry  of  247  ft.  Local  transmitter 
tests  were  not  performed  at  this  site,  but  later  analysis  of  the  data  in  the  laboratory  revealed  that 
the  transducer  resonance  was  about  336  Hz  at  receiver  signal  levels.  This  was  a  distinct  mismatch 
with  the  transmitter  frequency  of  308  Hz  when  the  narrow  band  character  of  the  devices  is  consid¬ 
ered,  and  this  points  up  a  serious  drawback  in  attempting  to  use  these  transducers  for  both  transmit¬ 
ting  and  receiving.  There  could  hardly  be  any  doubt  that  the  signal-to-noise  ratio  was  diminished  by 
perhaps  10  to  15  dB  due  to  this  detuned  condition. 


FIGURE  23.  IRON  MINE -TRANSMITTER  INSTALLATION,  2475  LEVEL 


A  wire  line  for  remote  triggering  was  dropped  down  a  vertical  vent  raise  (about  6  to  S  ft 
in  diameter)  which  joined  the  two  Ievelsat  a  point  some  1 50  to  250  ft  from  each  transducer  instal¬ 
lation.  This  did  provide  a  surface  wave  propagation  path  along  the  drift  and  vent  raise  walls  of  about 
600  ft  with  two  90-deg  turns.  However,  from  the  data  analysis,  it  is  believed  that  lirtle  energy  of 
significance  reached  the  receiver  by  this  path.  At  the  time  of  the  test,  unfortunately,  mine  machinery 
was  operating  within  about  200  ft  of  each  transducer  site  so  that  there  was  considerable  seismic  and 
acoustic  noise.  Nevertheless,  the  transmitter  was  triggered  a  number  of  times  and  signals  from  the 
RADC  transducer  recorded.  Triggering  was  very  erratic,  causing  it  to  be  surmised  that  the  wire  line 
had  been  damaged  by  mine  machinery  during  the  preceding  night  shift, operations.  Four  typical 
recorded  pulses  arc  presented  in  Figure  24,  the  latter  being  obtained  during  a  relatively  quiet  period. 
The  upper  trace  >n  each  case  is  the  trigger  and  the  lower  the  raw  received  signal.  The  transmit  pulse 
was  about  0.14  sec  long,  and  this  seismic  input  signal  was  assumed  to  be  nearly  identical  to  that 
obtained  in  the  limestone  quaTy  and  given  in  Figure  20.  Even  with  the  extended  tail  on  the  received 
pulse,  due  to  reflections,  refractions,  etc.,  a  data  rate  of  2  to  3  pulses/sec  appears  quite  practice. 

The  quiet  pulse  is  repeated  in  Figure  25.  stretched  out  timewise  as  an  aid  in  determining 
the  propagation  velocity.  Using  the  same  method  as  employed  for  the  quarry  data  and  assuming  the 
trigger  actually  occurred  on  the  first  transient  on  the  trigger  line,  the  propagation  velocity  was  found 
to  be  over  50.000  ft/sec,  an  obvious  absurdity  since  the  compression  wave  velocity  in  the  porphyry' 
as  calculated  from  physical  constants  o»  the  material  obtained  from  the  mine  people  indicated  an 
expected  range  of  19.000  to  23,000  ft/sec.  Accordingly,  the  envelopes  of  the  assumed  input  and  out¬ 
put  pulses  were  carefully  plotted  and  compared  with  theoretical  response  curves,  and  it  v.as  found 
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that  the  response  more  nearly  corresponds  to  that  of  two  first-order  elements,  with  time  constants 
differing  by  a  factor  of  four.  Thus,  the  time  required  for  the  input  to  reach  63  percent  of  final  value 
corresponds  to  the  time  required  for  the  output  to  reach  50  percent  of  final  value.  Using  this  crite¬ 
rion.  the  propagation  velocity  was  found  to  be  about  17.000  ft/sec.  This  agrees  closely  with  the  cal¬ 
culated  values  such  this*  propagation  by  compressional  waves  is  assumed.  (Note  that,  since  the  prop¬ 
agation  velocity  in  the  magnetite  ore  body  was  less  than  that  in  the  porphyry,  there  was  no  tendency 
to  propagate  a  guided  wave  up  this  interface:  see  Owen/1’) 

With  a  signal-to-noise  ratio  approaching  the  order  of  20  dB  demonstrated  (during  quiet 
periods  at  least),  it  was  decided  to  extend  the  range  and  attempt  further  communications  experi¬ 
ments.  These  will  be  described  next. 

c.  Receiver  Installation-  1825  Level 

A  receiving  site  on  the  1 825  level  was  selected.  This  site  was  642  ft  vertically  above  and 
406  ft  horizontally  offset  front  the  transmitter,  giving  a  slant  range  through  the  porphyry  of  760  ft. 
This  was  a  relatively  quiet  location  It  was  not  feasible  to  interconnect  the  iwo  sites  with  a  remote 
trigger  line:  therefore,  the  transmitter  was  pulsed  manually  by  prearranged  time  schedule.  Figure  26 
is  a  photograph  of  the  receiving  site  on  the  1825  level  showing  the  RADC  transducer  affixed  to  the 
rib  of  the  drift.  Again,  it  was  not  feasible  to  perform  local  transmitter  tests  at  this  site,  although 
later  laboratory'  analysis  of  the  data  indicated  that  the  installed  resonance  of  the  RADC  transducer 
was  only  260  Hz,  a  very  "‘sour”  installation  indeed  and  very  badly  detuned  with  respect  to  the 


transmission  frequency.  It  was  remarkable  that  signals  were  received  at  all  from  this  unit.  The  sen¬ 
sitive  accelerometer  was  also  used  at  this  site,  most  of  the  time  in  conjunction  with  a  filter  having  a 
100-Hz  pass  band  centered  on  330  Hz  with  24-dB/octave  skirts.  The  accelerometer  was  hand  heia 
against  the  rock  face  at  various  points  as  a  technique  to  maximize  the  received  signal  as  observed  on 
the  oscilloscope. 


Typical  signals  received  by  the  sensitive  accelerometer  and  the  RADC  transducer  are  shown 
in  Figure  27.  The  inputs  here  were  manually  gated  pulses  of  about  1 /2-sec  duration.  The  signal-to- 
noise  ratios  are  obviously  poor,  particularly  for  the  RADC  transducer,  although  most  of  the  raw  data 


FIGURE  27.  IRON  MINE-TYPICAL  RECEIVED  SIGNALS,  1825  LEVEL 

exhibited  even  poorer  signal-to-noise  ratios  as  will  be  shown  later.  It  is  of  considerable  interest  to 
note  that  the  gated  pulse  signal  was  audibly  detectable  by  the  unaided  ear,  though  at  very  low  level, 
throughout  most  of  the  test  period.  By  use  of  a  microphone,  a  short  sample  of  this  was  recorded  on 
tape.  After  2  days  of  tests  on  this  level,  the  raw  data  did  not  appear  to  j  jstify  a  further  extension  of 
the  transmission  range  (say  by  moving  the  receiver  site  to  a  still  higher  le/el  in  the  mine),  and  so  the 
experiments  were  terminated. 

3.  LABORATORY  SIGNAL  PROCESSING 

In  the  laboratory,  the  tape-recorded  signals  were  translated  through  a  32: 1  time  expansion  and 
played  out  on  a  chart  recorder  for  detailed  analysis  and  evaluation.  In  addition,  it  was  decided  to 
attempt  some  signal  processing  to  enhance  the  received  signal-to-noise  ratios.  Various  pass  band  fil¬ 
ters  were  tried  with  moderate  degrees  of  success,  until  the  data  were  played  through  the  narrow  pass 
band  filter  of  a  Hewlett  Packard  HP  302  A  wave  analyzer.  This  filter  has  a  pass  band  of  7  Hz  (—3  dB) 
and  70  Hz  (-60  dB).  The  results  achieved  were  rather  impressive,  and  typical  examples  are  presented 
in  Figures  28,  29,  and  30.  In  Figure  28,  the  RADC  transducer  signal  on  the  2275  level  is  displayed 
along  with  the  HP  302  A  filtered  signal.  The  improvement  in  signal-to-noise  ratio  is  remarkable  and 
approaches  40  dB  for  the  pulse  received  during  the  quieter  period.  Similarly,  Figure  29  presents 


FIGURE  29.  IRON  MINE-NARROW  BAND  FILTERING,  1825  LEVEL 


!i 


Ilillli 


illiil 


i^Mi 


FILTERED  SIGNAL 


i  i  m  it  i 

5/32  SEC  |-* 


^ 1 

gat 


_ _ _ hmhp" 

'inag»fS?Wr"F^i¥s'SW»TT^iFm’ 


liili 


fetssitaUtil 


RAW  SIGNAL,  SENSITIVE  ACCELEROMETER 


FIGURE  30.  IRON  MINE-NARROW  BAND  FILTERING,  1825  LEVEL 

typical  results  for  the  filtered  RADC  transducer  on  the  1 825  level.  The  signal-to-noise  ratio  in  the 
raw  signal  was  .ery  poor,  and,  again,  the  improvement  in  the  signal-to-noise  ratio  of  the  filtered  sig¬ 
nal  is  remarkable,  perhaps  exceeding  20  dB.  The  manually  gated  transmitter  pulses  were  about 
3/4  sec  in  duration  in  this  example.  The  effect  of  narrow  band  filtering  on  the  sensitive  acceler¬ 
ometer  signal  recorded  on  the  1825  level  is  shown  in  Figure  30.  The  record  shows  two  short  pulses 
(about  1/4  sec  each)  and  one  long  pulse  (about  1  sec),  such  as  might  be  used  in  a  coded  message.  It 
is  seen  that  the  narrow  band  filter  does  give  a  lagged  response,  but  the  rise  time  of  the  otherwise 
Q-limited  pulse  is  not  greatly  increased.  A  data  rate  of  two  pulses/sec  appears  to  be  quite  practical 
even  with  a  narrow  band  filter  of  the  type  employed  here. 

4.  SIGNAL  ATTENUATION  WITH  RANGE 


Relative  received  signal  levels  on  the  2275  and  1825  levels  were  compared  for  purposes  of  esti¬ 
mating  signal  attenuation  loss  in  the  seismic  medium.  The  loss  was  estimated  to  be  about  15  dB, 
although  the  accuracy  of  the  data  is  open  to  serious  doubt  since  the  RADC  receiving  transducer  was 
detuned  at  both  receiving  sites  to  differing  degrees.  Unfortunately,  because  of  equipment  difficulties 
at  the  time,  the  sensitive  accelerometer  signal  v/as  not  recorded  on  the  2275  level. 


. 


SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


Results  of  communications  experiments  in  an  operating  iron  mine  have  shown  that  the  RADC 
piezoelectric  transducer  performed  quite  well  as  a  source  or  transmitter  for  radiating  communica- 
tions-type  gated  pulse  signals  into  hard,  competent  rock  media.  For  a  peak  radiated  power  of  about 
30  W,  signals  were  received  at  a  distance  of  760  ft  through  hard  rock,  which,  after  being  processed 
through  a  narrow  band  filter  tuned  to  the  transmitter  unit  frequency,  exhibited  a  signal-to-noise 
ratio  of  the  order  of  20  dB  and  a  data  rate  capability  of  about  two  pulses/sec.  These  results  taken 
together  v/ith  the  results  of  the  theoretical  analysis,  laboratory  tests,  and  other  field  tests  lead  to  the 
following  conclusions: 

(1)  A  piezoelectric  transducer,  loaded  by  a  high  Q  spring-mass  resonator  on  the  air  side,  by 
proper  design  can  be  efficiently  coupled  (or  impedance  matched)  to  hard,  competent 
rock  media  for  the  purpose  of  radiating  seismic  signals  into  the  media. 

(2)  The  power  conversion  efficiency  of  such  a  source  transducer  can  be  made  about  70  to 
80  percent,  while  overall  power  conversion  efficiency  employing  a  class  B,  solid-state 
driver  can  be  made  to  approach  50  percent. 

(3)  Because  of  the  narrow  operating  bandwidth  of  such  a  liigh  Q  transducer,  pulse  width  or 
pulse  code  modulation  techniques  are  indicated,  and  data  rates  will  be  limited  to  a  few 
pulses  per  second. 

(4)  Transducer  resonant  frequency  of  about  300  Hz  is  probably  a  good  compromise  between 
power  radiating  capability  and  attenuation  loss  in  hard  rock  seismic  media.  However, 
definitive  data  to  support  this  statement  were  not  generated  in  this  effort. 

(5)  Such  a  transducer  it  not  suitable  for  use  as  a  receiver  because  of  the  variations  in  installed 
resonant  frequency  at  each  site.  A  sensitive  receiver  such  as  an  accelerometer  or  geophone, 
exhibiting  say  a  50-Hz  bandwidth,  would  be  a  better  selection. 

(6)  To  maximize  the  signal-to-noise  ratio,  particularly  in  noisy  environments,  the  received 
signalslould  be  filtered  by  a  very  narrow  band  filter  tuned  to  the  transmitter  unit  resonant 
frequency. 

(7)  Seismic  propagation  through  highly  fractured  media  will  result  in  large  attenuation  losses 
and  should  be  avoided  if  possible. 

(8)  Attachment  of  the  transducer  to  the  rock  medium  interface  should  be  done  very  carefully 
to  insure  good  coupling.  The  coupling  should  be  verified  by  performing  power  radiation 
tests  and  the  transducer  moved  and  reinstalled  if  indicated. 

(9)  In  massive  hard  rock,  effective  seismic  communication  ranges  of  1000  ft  or  greater  should 
be  practical  with  total  peak  power  consumption  of  100  W. 
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(10)  For  shorter  seismic  transmission  paths  tlirough  rock,  higher  frequencies  could  be  employed 
with  correspondingly  increased  data  rate  capabilities.  Depending  upon  the  frequency 
selected  and  the  signal  attenuation  in  the  rock,  this  might  require  more  or  less  power  radi¬ 
ating  capability  than  the  current  RADC  transducer.  For  sufficiently  short,  low-loss  paths, 
a  nonresonant,  broadband  piezoelectric  transducer  might  become  feasible. 

(11)  Adaptation  of  existing  technology  in  piezoelectric  transducers  and  in  electronic  techniques 
and  equipment  to  this  application  is  straightforward. 

The  above  conclusions  are  meant  to  demonstrate  the  feasibility  of  employing  this  type  of  piezo¬ 
electric  transducer  to  seismic  communications,  at  least  within  the  limits  specified.  Since  the  relative 
success  of  such  an  application  is  highly  dependent  upon  the  seismic  propagation  characteristics  of  a 
given  transmission  path  of  interest,  it  is  recommended  that  preliminary  communications  tests  be 
performed  at  any  such  site  before  finalizing  of  design  and  equipment  parameters. 

Further,  such  tests  could  be  greatly  facilitated  and  the  design  of  a  final  transducer  further 
advanced  by  developing  an  improved  version  of  the  current  RADC  transducer.  Specifically,  the 
objectives  for  such  a  development  would  be  as  follows: 

( 1 )  Increased  radiating  power  capacity ; 

(2)  Elimination  of  nonlinearity  problem; 

(3)  Reduction  in  size  and  weight: 

(4)  Design  of  a  ruggedized  electronics  drive  and  receiving  system. 


APPENDIX 


DERIVATION  OF  PIEZOELECTRIC  CRYSTAL  EQUIVALENT  CIRCUIT 


Consider  a  piezoelectric  crystal  of  thickness,  £,  and  cross  sectional  area,  a,  driven  in  the  thick¬ 
ness  mode,  and  neglect  radial  and  other  modes.  As  given  by  Gooberman(l8)  or  Mason(l9),  the  appli¬ 
cable  piezoelectric,  equations  are 


Tx-Cx  ~~  ~  hDx 
ox 


(22) 


Ex 


(23) 


where 

Tx  —  tension  in  crystal  along  the  x  direction,  N/m2 

Cx  —  elastic  constant  along  the  x  direction  at  constant  charge  density,  N/m2 
£*  -  crystal  displacement  along  *  direction,  m 

h  —  piezoelectric  constant,  V/m  or  N/C 
Dx  -  charge  density  in  x  direction,  C/m2 
Ex  -  electric  field  in  x  direction,  V/m 

ei  -  absolute  permittivity  in  x  direction  at  constant  material  strain,  F/m 

MKS  units  are  shown  since  this  system  of  units  is  somewhat  more  convenient  to  work  with. 

Eqs.  (22)  and  (23)  express  the  relationships  among  the  electrical  and  mechanical  conditions  at 
a  given  point  within  the  crystal.  If  the  propagational  wavelength  is  very  long  compared  to  the  crystal 
thickness,  £  along  the  x  direction  (such  as  is  the  case  here;  i.e.,  the  wavelength  in  PZT-4  material  at 
350  Hz  is  about  38  ft),  then  the  equations  become,  upon  dropping  the  subscript,  x, 

T  =  CP  —  —  hD 
£ 

1  A£ 

E=-D~h — 
e*  £ 


where  A£  is  the  total  expansion  in  crystal  thickness.  Multiplying  the  first  by  the  area,  a,  gives 


F=  kAZ—hQ 


(24) 


where 


CDa 

k  =  -  — ,  the  short  circuit  spring  rate  of  the  ciystal 


Q  =  Da,  free  charge  on  ciystal 
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and  F  is  the  tension  force  across  ciystal.  Multiplying  the  equation  for  £  by  the  thickness,  2  gives 


V  =  ^--hAS. 


where  V  is  the  voltage  across  crystal  and 


C0  =  — ,  clamped  capacitance  of  crystal 
Combining  Eqs.  (24)  and  (25)  to  eliminate  Q  and  using  a  compressive  force 


yields 


Fr=-F 


Fc=hC0V-At(k-h2C0) 


For  various  reasons,  it  is  convenient  and  customary  to  work  with  the  concept  of  impedance, 
both  electrically  and  mechanically,  where  electrical  impedance  is  the  ratio  of  voltage  to  current  and 
mechanical  impedance  is  the  ratio  of  force  to  velocity.  This  implies  the  following  analogical  cor¬ 
respondence: 


and,  derivatively, 


voltage  — *  force 
current  — *•  velocity 
charge  — *•  displacement 


capacitance  — *■  inverse  spring  rate 
inductance  — ►  mass 
resistance  7— ►  viscous  damping 


Thus,  if  sinusoidal  time  variation  at  radian  frequency,  co,  is  assumed  and  velocity  is  defined 


and  current  is 


v=jcoAi 


i=juQ 


then  Eqs.  (25)  and  (26)  become,  upon  some  rearranging, 

i  =  joiCQ  V  +  hC0  v 


Fc  =  hCQ  V  —  —  (k  —  h2Ca ) 
/w 


These  equations  suggest  the  following  equivalent  circuit: 
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V.hC0 
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Ideal  Electrical-Mechanical 
Transformer 
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All  quantities  to  the  left  of  the  transformer  are  electrical,  while  all  quantities  to  the  right  are  mechan¬ 
ical.  The  mechanical  “capacitors”  are  proportional  to  inverse  spring  rates,  one  of  which  represents 
the  elastic  spring  rate  of  the  crystal  when  electrically  shorted  and  the  other  a  negative  spring  rate 
associated  with  the  change  in  free  charge  resulting  with  constant  voltage  excitation.  These  two 
capacitors  shall  be  combined  into  a  single  capacitor  in  what  .ollows. 

For  n  crystals,  stacked  mechanically  in  a  series  but  connected  electrically  in  parallel,  the  follow¬ 
ing  equivalent  circuit  is  suggested,  where 


Cc  k-h2Ca 
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It  is  evident  that  the  total  current  is 


i?  =  ni 


the  output  velocity  is 


Vf  =  nv 


and  the  output  force,  Fc,  across  the  entire  stack  is  the  same  as  the  force  across  a  single  crystal.  Con¬ 
sequently,  Eqs.  (27)  and  (28)  become,  for  a  stack  of  n  crystals, 

iT  =  junC0  V  +  hC„vT  (29) 

Fc  =  hC0  V  — —(k  —  h2C0)  (30) 

icon 

and  thr  fxnai  equivalent  circuit  is  given  as  follows  where  we  have  now  added  a  mechanical  load  to 
each  end  of  the  crystal  stack  represented  as  mechanical  impedances  Z,  and  Z2  and  have  dropped  the 
subscripts  on  v  and  i: 


n 


Alternatively,  the  mechanical  portion  of  the  circuit  might  be  depicted  by  using  mechanical 
elements  as  shown  on  the  following  page. 

From  manufacturers’  published  data,  numerical  constants  for  the  crystal  stack  are  as  follows: 

n  =  6 

a-  2.85  in.2  =0.00184  m2 
8  =  0.2  in.  =  0.00508  m 
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